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Abstract. We obtain several results on the computation of different and dis- 
criminant ideals of finite extensions of local fields. As an application, we deduce 
routines to compute the p-adic valuation of the discriminant Disc(/), and the 
resultant Res(/, g), for polynomials f(x),g(x) £ A[x], where A is a Dedckind 
domain and p is a non-zero prime ideal of A with finite residue field. These rou- 
tines do not require the computation of neither Disc(/) nor Res(/, g); hence, 
they are useful in cases where this latter computation is inefficient because the 
polynomials have a large degree or very large coefficients. 



Introduction 

Let A be a Dedekind domain whose field of fractions K is a global field. Typical 
instances of A are the ring of integers of a number field, or the polynomial ring in 
one indeterminate over a finite field. Let L/K be a finite separable extension and 
B the integral closure of A in L. Let 9 G L be a primitive element of L/K, with 
minimal polynomial f(x) S A[x]. Also, let p be a non-zero prime ideal of A, Vp the 
canonical p-adic valuation and Kp the completion of if at p. 

The Montes algorithm 5 computes an OM representation of every prime ideal 
*P of B lying over p [7] . An OM representation is a computational object supporting 
several data and operators, linked to one of the irreducible factors (say) F{x) of f(x) 
in Kp [x] . Among these data, the OM representation contains the Okutsu invariants 
of F, which reveal a lot of arithmetic information about the finite extension of Kp 
determined by F [TT1I4]. 

In this paper, we present an algorithm to compute the ^P-adic valuation of the 
different ideal Diff(L/iv') in terms of the OM representation of *p. Of course, this 
determines the p-adic valuation of the discriminant Disc(L / K) . This computation 
does not require the computation of the discriminant Disc(/) of f(x); actually, 
the algorithm provides the value of u p (Disc(/)) as a by-product. We use this fact 
to derive a routine to compute u p (Disc(g)) for an arbitrary polynomial g(x) £ 
K[x]. Also, we use similar ideas to compute the p-adic valuation of the resultant 
Kes(g,h) of two polynomials g(x),h(x) £ K[x]. These routines may be useful in 
cases where the natural computation of Disc(g) or Res(g, h) is inefficient because 
the polynomials have a large degree or very large coefficients. 

The algorithms are described in sections 2, 3. If A/p is small, the compu- 
tation of 5 := v p (Disc(g)) requires O (n 2+e + n 1+e S 2+e ^j word operations, where 
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n = dcgg. The computation of v p (Res(g, h)) has the same complexity, but tak- 
ing n — max{degg,degft.}, S — v p (Res(g, h)). The algorithms have an excellent 
practical performance. In section 4 we present some numerical tests with Magma, 
in which we compare running times with the naive algorithms that compute first 
Disc(<?) or Res(g, h), and then its p-valuation. 

The core of the paper is section 1, where we discuss the computation of the 
different ideal of an extension of local fields. Let L<$ be the completion of L at the 
prime ideal *p. Consider an embedding, L C £<p <-> K p , of L into a fixed algebraic 
closure of K p , and let F(x) G ^"p[^] be the monic irreducible factor of f(x) that 
vanishes on the image of 9 in K p , which we denote still by the same symbol 9. 
Let e(F) = e(ty/p) be the ramification index of *P over p, and f(F) — /(5}3/p) the 
residual degree. The different ideal Diff(Ltp/ K p ) is equal to the (e(F) — 1 + p)-th 
power of the maximal ideal of the ring of integers of Ltp , where p is a non- negative 
integer that vanishes if and only if is tamely ramified over p. On the other hand, 
from Theorem 11.71 we deduce identities: 

Up (Disc(F)) = nv p (F'(d)) = np(F) + f(F)p, 

for some invariant p(F) G Q. The three numbers e(F), f(F), fi(F) are Okutsu 
invariants of F, which can be expressed in terms of data contained in the OM 
representation of J 7, by a direct formula (Proposition II .4[) . Therefore, for the com- 
putation of Diff (Lrp/K p ) and u p (Disc(i 7 ')), we need only to compute v p (F'(9)). 
Since it is impossible in general to get an exact computation of the polynomial 
F(x), in Theorem 11.111 we show how to construct a polynomial 4>(x) G A[x] which 
is sufficiently close to F(x) to have v p ((f>'(9)) = v p (F'(9)). Finally, once we get this 
sufficiently good approximation <fi to F, the value of v p ((f>'(9)) may be determined 
by using Newton polygons of higher order (Proposition 1 1 . 5[) , whose computation 
relies as well in data contained in the OM representation. 

1. Computation of local different ideals 

Let k be a local field, i.e. a locally compact and complete field with respect to a 
discrete valuation v. Let O be the ring of integers of k, m the maximal ideal, and 
F = O/m the residue field, which is a finite field. Let p be the characteristic of F. 

Let fe scp C k be the separable closure of k inside a fixed algebraic closure. Let 
v. k — > Q U {oo}, be the canonical extension of the discrete valuation v to k, 
normalized by v(k) — Z. 

1.1. Okutsu invariants of an irreducible separable polynomial. Let F(x) £ 
0[x] be a monic irreducible separable polynomial, 9 G fc scp a root of F(x), and 
L = k(9) the finite separable extension of k generated by 9. Denote n :— [L: k] — 
deg F. Let Ol be the ring of integers of L, mj the maximal ideal and F^ the residue 
field. We indicate with a bar, : 0[x] — > F[x], the canonical homomorphism of 
reduction of polynomials modulo m. 

Let ... , 4> r ] be an Okutsu frame of F(x), and let </V+i be an Okutsu approxi- 
mation to F(x). That is, </>i, . . . , r +i G 0[x] are monic separable polynomials of 
strictly increasing degree: 

1 < mi :— deg 0i < • ■ • < m r :— deg0 r < m r+ i :— deg0 r+ i = n, 
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and for any monic polynomial g(x) G 0[x] we have, for all < i < r: 

^ v(g(9)) v(M0)) < v(A+i(e)) 

mi < deg.g < m i+1 => — < < , 

degg m l m i+1 

with the convention that mo = 1 and (j)o(x) — 1. It is easy to deduce from these 
conditions that the polynomials <fii(x), . . . , (f> r+ i(x) are all irreducible in 0[x\. 

The length r of the frame is called the Okutsu depth of F(x). Okutsu frames were 
introduced by K. Okutsu in [11] as a tool to construct integral bases. Okutsu ap- 
proximations were introduced in [3], where it is shown that the family fix, ... , <j) r +i 
determines an optimal F -complete type of order r + 1: 

(1) t F = (V>o;(<£i,Ai,V>i);--' ;(0 r ,A r ,'0 r );(0 r+ i,A r+ i,-0r+i))- 

In the special case 4> r +i = F, we have A r +i = — oo and ?pr+i is not defined. We 
call tp an OM representation of F. 

An OM representation of the polynomial F carries (stores) several invariants and 
operators yielding strong arithmetic information about F and the extension L/k. 
Let us recall some of these invariants and operators. 

Attached to the type tp, there is a family of discrete valuations of the rational 
function field k(x), the MacLane valuations: 

v,: k(x) — >ZU{oo}, l<i<r + l, 

such that = vi(F) < ■ ■ ■ < v r+ i(F). The wi-value of a polynomial in k[x] is the 
minimum of the u-values of its coefficients. 

Also, tp determines a family of Newton polygon operators: 

Ni:k[x] — >2 R \ l<i<r+l, 

where 2 R is the set of subsets of the Euclidean plane. Any non-zero polynomial 
g(x) G k[x] has a canonical ^-development: 

g(x) = a s {x)4>i(x) s , dega s <m l; 

A — '0<s 

and the polygon Ni(g) is the lower convex hull of the set of points (s,Uj(a s 0|)). 
Usually, we are only interested in the principal polygon N~ '(g) C Ni(g) formed 
by the sides of negative slope. For all 1 < i < r, the Newton polygons Ni(F) 
and Ni((j)i + i) are one-sided and they have the same slope, which is a negative 
rational number \ e Q<o- The Newton polygon N r +i(F) is one-sided and it has 
an (extended) integer negative slope, which we denote by A r+ i E Z <0 U {— oo}. 

There is a chain of finite extensions: F = Fo C Fi C • • • C F r+ i = F^. The 
type tp stores monic irreducible polynomials ipi(y) £ ^i[y] such that Fj+i ~ 
¥i[y]/(ipi(y)). We have tpi(y) ^ y, for all i > 0. Finally, there are also residual 
polynomial operators: 

R l : k[x] — >Wi[y], 0<i<r + l. 

For all < i < r, we have Ri(F) ~ and Ri((pi + i) ~ fy, where the symbol ~ 

indicates that the polynomials coincide up to a multiplicative constant in F*. For 
% = we have Ro(F) = F = -0O 1 an< i Ro(<l>i) = ffri = ipo- The exponents uj i+ i are 
all positive and tu r +i = 1- The operator R r +i is defined only when 4> r +i ^ F; in 
this case, we also have R r +i(F) ~ Vv+i, with ifj r +i(y) G F r +i[y] monic of degree 
one such that ip r +i(y) ^ y. 
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We recall that the length of a Newton polygon N is the abscissa of its right end 
point; we denote it by £{N). There is a strong link between the operators Ri and 
N~ +1 Lem. 2.17,(2)]. 

Lemma 1.1. Let < i < r. For any non-zero polynomial g(x) G k[x], (-{N~ +1 {g)) 
is equal to the maximal exponent ord^, i Ri(g) with which ipi divides Ri{g) in Fj[y]. 

From all these data of the OM representation some more numerical invariants 
are deduced. Initially we take: 

to := 1, / :=deg?/'o, e := 1, h := V := [i := is = Q. 

Then we define for all 1 < i < r + 1: 

hi, ei positive coprime integers such that A^ = —hi/ei, 

f t . := degipi, 

m !: := deg<^ = ej_i/i_imi_i = (e e x ■ ■ ■ ei_i)(/ /i • • • fi-i), 

l l i ■= 111 r . :< Jj ■ • ' e ifi ~ !) V( e l • ■ ■ e j): 
v i : =El<i<i^'/( e l--- e i). 

Vi := = e^ifi-ifa-iVi-i + h^x) = (e • • • ei_i)(^_i + ^-i). 

The general definition of a type may be found in [51 Sec. 2.1]. In later sections, 
we shall consider types which are not necessarily optimal nor F-complete. So, it 
may be convenient to distinguish these two properties among all features of a type 
that we have just mentioned. 

Definition 1.2. Let t = (^o; {<j>i, Ai,^i); • • • ; {4>i, Xi,ipi)) be a type of order i. 

• We say that t is optimal if mi < ■ ■ ■ < rtii. 

• We say that t divides the polynomial g{x) € K[x] (and we write t | g) ifipi \ Ri{g) 
in Fj[y]. We say that t is g-complete i/ord^ ; Ri{g) = 1. 

Thus, for a general type of order i we have mi | • • • | m t and uii > 0, but not 
necessarily toi < • • • < m; = degF, and Ui = 1. These were particular properties 
of our optimal and incomplete type tp of order i = r + 1, constructed from an 
Okutsu frame and an Okutsu approximation to F. 

An irreducible polynomial F admits infinitely many different OM representa- 
tions. However, the numerical invariants et, f%,hi, for < i < r, and the MacLane 
valuations v x , . . . , v r +i attached to t^, are canonical invariants of F. 

The data A r +i,?/v+i are not invariants of F; they depend on the choice of the 
Okutsu approximation 4> r +i. The integer slope A r +i = —h r+ i measures how close 
is 4> r +i to F . We have 4> r +\ = F if and only if h r+ \ = oo. 

Definition 1.3. An Okutsu invariant of F(x) is a rational number that depends 
only on ei, . . . , e r , /i, . . . , /,-, hi,..., h r . 

Several arithmetic invariants of the polynomial F(x) and the field extension L/k 
are Okutsu invariants of F(x). For instance, let us define: 

e(F) :— e(L/k), the ramification index of L/k, 
f(F) := f(L/k), the residual degree of L/k, 

:= ma,x{v(g(9)) \ g(x) £ 0[x] monic of degree less than n}, 
ind(F) := length (0 L /0[6l]), (length as an O-module), 
f(F) := mm{6 € Z> | m s L L C O[0}}, 
exp(F) := min{<5 6 Z> | m s O L C G[6]}. 
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We abuse of language and say that ind(F) is the index of F(x), f(F) is the con- 
ductor of F(x), and exp(F) is the exponent of F(x). Actually, f(F) is the triL-adic 
exponent of the conductor of 0[8] as an order of Ol, and exp(F) = \f(F)/e(F)~\. 
Also, if k is a finite extension of the field Q p of p-adic numbers, then ind(i^) is re- 
lated to the u-adic exponent of the finite index (Ol ■ O[0\) of these abelian groups; 
more precisely, ind(F) = v((O l : 0[6)))/[k: Q p ]. 

The next proposition gathers the results [6l Cor. 3.8], Thm. 5.2] and 
Prop. 3.5]. 

Proposition 1.4. 

e(F) = e ei • ■ -e r , 
f(F) = foh---f r , 

= fl r = El<i<r( e i/7 ' ' ' e rfr ~ (ei • • • ej), 

exp(F) = [Mr J 

md(F) = |(/x r -l + e (F)- 1 ). 

We end this section with a result extracted from [51 Prop. 3.5,(5)]. 

Proposition 1.5. Let g(x) € 0[x] a non-zero polynomial. For any 1 < i < r take 
a line of slope Xi far below Ni(g), and shift it upwards till it touches the polygon 
for the first time. Let (0, H) be the intersection point of this line with the vertical 
axis. Then, 

v(g(0)) > v i+1 {g)/(e 1 ■■•ei) = iT/(e ei • • -e 4 _i), 
and equality holds if and only if ipi \ Ri(g) in Fj[y]. 

We emphasize two particular instances of this result, that occur quite often. 

Corollary 1.6. With the above notation, 

(1) Ifdegg < m i+ i, then v(g{9)) = v i+ i(g)/(e x ■ ■ ■ e^), for all 1 < i <r. 

(2) v{(j>i{e)) = (Vi + \Xi\)/{e ei ■ • ■ e^i) = w _! + v u for all 1 < i <r + I. 

1.2. Different and discriminant. Let Diff(L/fc) be the different ideal of L/k 
and Disc(L/fe) = ^/^.(Diff (L/k)), the discriminant ideal of L/k. From now on we 
denote 

e := e (F) = e(L/k), f := f(F) - f(L/k). 
The following identities are well-known: 

Disc(F) - Res(F, F') = N L/k (F'(9)), 
(2) v{Bisc(F)) = 2 ind(F) + u m (Disc(L/fc)), 

Diff(i/fc) = (m L ) e ~ 1+p , p > 0, 

and p = if and only if L/k is tamely ramified. 

Since the Okutsu invariants of F can be computed from data contained in the OM 
representation of F, the only extra information we need to compute the different 
ideal is the value of the non-negative integer p, which is a measure of the wildness of 
the ramification. If L/k is a Galois extension, then it is well-known how to deduce 
the value of p from the higher ramification subgroups of the Galois group [15] . We 
are interested in the computation of p for an arbitrary extension L/k. 

By applying the formula for ind(i 7 ') in Proposition [L4l we get a first formula for 
p in terms of the Okutsu invariants p(F), f(F) and the v- value of Disc(F). 



6 



NART 



Theorem 1.7. v(Disc(F)) = f(efi(F) + p) = np{F) + fp. 

However, we are interested in a computation of p that does not require the 
computation of Disc(i 7 '). This can be achieved as follows. 

Corollary 1.8. v r+1 (F') = v mL (F'(9)) = ep{F) + p. 

Proof. Since degF' < n = m r+1 , item 1 of Corollary 11.61 and Proposition 11.41 show 
that v(F'(9)) = v r+1 (F')/e. Hence, by © and TheoremO 

nv{F'{9)) = v(Disc(F)) = np{F) + fp. 

If we divide both terms of the equality by /, we get the result. □ 

As a consequence, the conductor of F(x) is an Okutsu invariant too. In fact, it 
is well-known that QUI Ch.III,§6,Cor.2] 

F'(6)O l = (m L )« F ) Diff(L/fc). 

On the other hand, Corollary 11.81 shows that 

v mL (F'(9)) = en(F) +p = ep(F) + v mL (Diff (L/k)) -e + l. 

Thus, by comparing the two identities we get an expression of f(F) as an Okutsu 
invariant. 

Corollary 1.9. f(F) = ep(F) -e + l = 2ind(F)//. 

In practice, we deal with some polynomial f(x) with coefficients in a Dedekind 
domain A, and F(x) is an irreducible factor of f(x) over the completion of A at 
a nonzero prime ideal p of A. Thus, we never get an exact computation of the 
polynomial F(x), but only an approximation to it. Therefore, if L/k is wildly 
ramified and we want to apply Corollary 11.81 to compute p, we need to find a 
polynomial (f>(x) £ A[x], sufficiently close to F(x) to have v r+ i(4>') = v r+ i(F') and 
p(4>) = p{F). To this end is devoted the next section. 

1.3. Okutsu approximations and the different. Suppose L/k is wildly rami- 
fied. Let (j)(x) £ 0[x\ be an Okutsu approximation to F{x). That is, is a monic 
separable polynomial of degree n such that v((f>(9))/n > v(4> r (9))/m r , where r is 
the Okutsu depth of F(x) and [fa,..., <fi r ] is some Okutsu frame of F(x). This 
polynomial <p(x) is always irreducible and it has the same Okutsu invariants as 
F(x) [U Lem. 4.3]. In particular, Proposition II .41 shows that 

(3) e(4>) = e(F), /(<£) = f{F), p{(j)) — p{F). 

However, the non-negative integer p(F) := v r+ i(F')~e(F)^(F) is not an Okutsu 
invariant of F(x), and we may have i> r +i(</>') ^ v r+ i(F'). The aim of this section is 
to find how close to F must be taken the approximation 0, to ensure that v r +i(4>') = 
v r+ i(F'). By Corollary 11.81 and ([3]), this implies p(<fi) = p(F); thus, the polynomial 
</> may be used to determine the value of p(F) and hence the different ideal of L/k. 

For any Okutsu approximation we may apply the results of section 11.11 to 
4> r +i := 4>. Since deg</> = deg-F and both polynomials are monic, the difference 
a(x) :— F(x) — 4>(x) has degree less than n = m r+ i. Hence, F{x) — a{x) + <p(x) is 
the canonical 0-development of F, and the right end point of N r+ i(F) is (1, V^+i). 

Since ord^, r R r (F) = w r+ i = 1, Lemma [TTTI shows that £(N~ +1 (F)) = 1. Hence, 
N~ +1 (F) = N r+ i(F) is one-sided of slope A r+ i = — h^, for some positive integer 
hcj, := w r +i(a) — V r +i (see Figure[I|. In particular, h r+ i — h^, e r +\ = 1, f r +\ = 1. 
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FIGURE 1. Newton polygon N r+1 (F) 



«r+l(o) 
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Clearly, v(<f>(9)) is a measure of the quality of <j> as an approximation to F. By 
item 2 of Corollary II .61 we have 

v(<t>(9)) = {V r+ i + h<j,)/e = /x r + i/ r +i- 

Since and e are constant, the positive integer is a measure of the quality 
of the approximation too. The main result in this section (Theorem II . lip finds 
a sufficient condition to ensure that <f> is sufficiently close to F to have the same 
exponent of the different ideal. Before proving this theorem we need a lemma, 
inspired by [P2] . 

Lemma 1.10. Let F(x) £ 0[x] be a monic irreducible separable polynomial, tp an 
OM representation of F as in (QJ), and 9 £ fc scp a root of F. 

Let g(x) £ 0[x] be a polynomial satisfying degg < m^+i, for some < i < r. 
Then, v(g'(9)) > v(g(9)) — z/j. 7/p | eo ei • • • ej_i, i/ien i/iis inequality is strict. 

Proof. Let us prove the lemma by induction on i. Suppose i — 0. Recall that 
= and _F = (ipo)" 1 , where ip is irreducible of degree / = mi. Let 7r £ m be a 
uniformizer, and express g(x) — -K a G(x), with G^O. Since degG < degg < mi = 
degV'Oj w e have ipo \ G, and this implies v(G(9)) = 0. Since g'(x) = ir a G'(x), we 
have then v(g'(9)) > a — v{g{9)). For i — 0, the empty product eo ei ■ • -e,_i is 
equal to one; hence, the condition p | eo e\ ■ ■ ■ ei-\ is never satisfied. 

Let i > 1, and suppose the statement of the lemma is true for all polynomials of 
degree less than m*. Assume mj < degg < m^i, and let 



be the canonical ^-development of g{x). 

For any s > 0, let (0, H s ) be the intersection point of the vertical axis with the 
line of slope passing through (s,Vi(a s <j)f)). Clearly, dega s (4>i) s < degg < mj+i; 
hence, Proposition II .51 and Corollary 11.61 show that 



and v(g(9)) = Vi+i(g)/e — H/(eo e\ ■ ■ ■ ej_i), where (0, H) is the intersection point 
of the vertical axis with the line of slope Xi that first touches Ni(g) from below (see 
Figure [2]). Therefore, 





v(a s (9)M8Y) = Hs/ieoe! ■ ■ ■ e,_i), 



(4) 



v(g(9)) 



e ei • • • e,_i 



H 



= mm 

0<s 



{ 



e ei • • • e l -\ 



} 



mm{v(a s (9)M0Y)}- 
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FIGURE 2. Computation of v i+ i 




Let us write: 

(j)i(x) g'(x) =y^ b s (x) 4>i(x) s , b s (x) = a' a (x) <f>i(x) + sa s (x) $ (x). 

* — »0<s<m 

Note that this is not the canonical ^-development of (j>i{x) g'(x). We claim that 

(5) v{b s (6)) > v(a s (6)) + Mi-i, for all s > 0, 

and if p eo e± ■ ■ ■ e»_i, then the inequality is strict. 

In fact, let us check that both summands of b s (6) = a' s (9)(f>i(9) + sa s {ff)4>' i {ff) 
satisfy this inequality. By the induction hypothesis and Corollarv ll. 61 (2). 

v{d s {9)4>i(6)) > v(a s (9)) - v^x + Hi-l + n > v(a s (6)) + m-i. 

On the other hand, [<j>\, . . . , 4>i-i\ is an Okutsu frame of (f>i(x) j4] Cor. 2.5], and 
Corollary 11.81 shows that Uj(<^) = eo e% ■ ■ • et-ipi-i + p((f)i). Since deg^ < m,, 
Corollarv ll. 61 shows that v((f>' i (9)) = + p(cf)i)/(eo e,\ • ■ ■ e-i-i)- Hence, 

v(sa s (9)(j)' i {9)) > v(a s (9)) + p t -i + p{<t>i)/{e ei ■ ■ ■ e^i) > v(a s (9)) + fM-i, 

and the inequality is strict if p | eo e\ ■ • ■ = e(4>i)- This ends the proof of (|5"T). 
Now, from J5]) and ((U), we get 

v{cj)i{e)g' '(9)) > min{v(b s (9)M9) s )} 

0<s 

> mm{v{a s {9)<p i (9) s )} + = v(g(9)) + Mi-i, 

0<s 

and the inequality is strict if p | eo e\ ■ ■ ■ e,_x- This ends the proof of the lemma 
because v(4>i(9)) = p,i-i + by Corollary II .61 □ 

Theorem 1.11. Let <j) be an Okutsu approximation to F, with > p(F). Then, 
p(4>) = p(F). In particular, if L' /k is the extension determined by (j), then Diff (L/k) 
and Diff(L'/fc) are divided by the same power o/rriL and m^', respectively. 

Proof. Let F(x) = a(x) + (f>(x), with a(x) G 0[x\ of degree less than n = m r+ i. 
Clearly, v(a(9)) = v(4>(9)) = p r + v r +i- Hence, Lemma \1 . 101 shows that 

v(a'(9)) > v(a(9)) - v r = p r + v r+x - v r = p r + — > p r + ^— ^ = v(F'(d)), 

e e 

the last equality by Corollarv ll .81 From F'(x) — a'(x)+4>'(x), we deduce v{4>' (9)) — 
v(F'(9)), or equivalently (by Corollary II .6[) , v r+ i((f>') = v r+ i(F'). 

Since (j) and F have the same Okutsu invariants, we have e{(f) — e(F) and 
p{cj)) = p(F). By Corollary Owe deduce that p{4>) = p{F). □ 
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In order to get an algorithm to compute p(F) (and Diff (L/k)), we must replace 
the condition h$ > p(F) by something checkable. This is done in a different way 
for local fields of equal or unequal characteristic. 

Corollary 1.12. If k has characteristic zero and > ev(e), then p(<j)) = p(F). 

Proof. It is well-known that in this case p(F) < ev(e). Hence, h$ > ev(e) implies 
h<t> > p(F), and the arguments of the proof of Theorem 11.111 show that v r +i{<fi') > 

V r+ l{F'). 

Let us now exchange the role of cf> and F. The list [4>i,- ■ ■ ,4> r ] is an Okutsu 
frame of <f>{x), and F(x) is an Okutsu approximation to 4>(x) [U Lem. 4.3]. Hence, 
v r+i(F) = V r+ \. Also, 4> = —a + F is the canonical F-development of (j>, so 
that hp — v r+ i(a) — V r+ \ = h<t, > ev(e). By the previous argument we have 
v r+i(F') > v r+ i((j)'), and hence, v r+ i(F') = u r +i (</>')■ ^ 

In the equal characteristic case, the value of p(F) is unbounded among all poly- 
nomials F of a fixed degree. We may proceed as follows in this case. 

Corollary 1.13. If k has characteristic p and > p(4>), then p((f>) = p{F). 

Proof. The argument of the proof of Theorem 11.111 shows that 

v(a'(9)) >p jr + ^± >f x r + B^l = „(0'(0)). 

e e 

Hence, v(<j)'(6)) = v(F'(0)), and this implies p(<f>) = p(F). □ 

Therefore, we may apply the Single-factor lifting (SFL) algorithm of [8] to im- 
prove the Okutsu approximation cf> till the checkable condition of either Corollary 
11.121 or Corollary 11.131 is satisfied, and then we obtain a computation of p(F) as 
p(4>). In each iteration, from a given Okutsu approximation 0, we compute a better 
one $ satisfying /i$ > 2/i^; thus, only a finite (and small) number of iterations are 
required. 

In the equal characteristic case, the condition h$ > p(4>) must be attained after 
a finite number of iterations, because p(F) is fixed and when we reach > p(F), 
then we have already p((j>) = p(F) < h^, by Theorem II. Ill 

2. Computation of the p-adic valuation of discriminants 

Let A be a Dedekind domain whose field of fractions K is a global field. Let p be 
a non-zero prime ideal of A, and v p : K — > ZU{oo} the canonical p-adic valuation. 
Let L — K[x]/ (f(x)) be the finite extension of K determined by a monic irreducible 
separable polynomial f(x) € A[x], and let B be the integral closure of A in L. 

In this section, we apply the results of the preceding section to compute the 
p-adic valuation of the discriminant ideal Disc(L/K), and the i> p -value of the dis- 
criminant Disc(g) of an arbitrary polynomial g{x) G Both tasks are based 
on a combination of the Montes algorithm [S] and the Single-factor lifting (SFL) 
algorithm [8]. 

2.1. Local computation of the different. Let K p be the completion of K at 
the prime ideal p and let O p be the ring of integers of K p . We denote still by v p 
the canonical extension of the p-adic valuation to an algebraic closure of K p . Let 

f(x)=F 1 (x)---F t (x), 
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be the factorization of f(x) into a product of monic irreducible polynomials in 
Op [a;]. Let Li, . . . ,Lt be the finite extensions of K p obtained by adjoining to K p 
a root of Fi,...,Ff, respectively. After a classical theorem of Hensel, there are 
exactly t prime ideals *Pi, . . . , *p t in B lying over p, and 

eflti/p) = e(L t /K p ), f^/p) = f(Li/K p ), l<i<t. 

Let us denote by Ftp (a;) £ Op[x] the irreducible factor that corresponds to any 
such prime ideal *}3 dividing p, and let Ltp be the corresponding extension of K p . 
At the input of f(x) and p, the Montes algorithm computes an OM representation 
of every *}3 | p , or equivalently, of every irreducible factor Ftp [7] : 

(6) t<p = (-0o; (0i, Ai, ip%); ■ ■ • ; (0 r , K, ipr)l (</v+i, A r +i, Vv+i)), 

carrying all invariants and operators described in section 11.11 The polynomials 
0i, ... , r +i have all coefficients in A. Each type ttp is /-complete and Ftp-complete; 
it singles out *}3 (or Ftp) by the following property: 

t v |Fp, ttpfFo, VQ^qj. 

In particular, all Okutsu invariants of Ftp, like e(Fp), /(Ftp), //(Fp), ind(Fp), 
f(Fp) and exp(Fp), are computed by direct formulas in terms of data supported 
by the OM representation of 

Remark 2.1. (1) Let 4>{x) £ A[x] be an Okutsu approximation to F<g(x), and 
let N r+ i be the Newton polygon operator with respect to 4>{x) and the valuation 
iV+i- By [6l Thm.3.1], we have N~ +1 {f) = A r+ i(Fp). Thus, the local invariant 
h^, that was defined by :— w r +i(Fp — 0) — V r +i, may be computed as well as 
h,), = v r +\(ao) — V r +\, where ao is the 0-th coefficient of the canonical 0-development 
of f(x) (see Figure [J). 

(2) The restriction to k of the valuation v r+ i coincides with e(^P/p)uqj. Hence, 
the value of w r +i(ao) does not change if we replace 0(x) by <f>(x) + b(x), with 
b(x) e A[x] of degree less than n«p and 6(x) = (mod p m ), for some m > (h$ + 
T0+i)/eCP/p). In particular, the coefficients of the Okutsu approximation may 
be always simplified modulo such a power of p. 

At the input of an OM representation of (or Ftp ) and a prescribed precision 
i/£N, the SFL algorithm computes an Okutsu approximation 0(x) to Fp(cc), such 
that <p(x) = F(x) (mod p v ) f8j. In each iteration of the main loop of SFL, a given 
Okutsu approximation <j){x) is used to construct a better Okutsu approximation 
$(x) with /i$ > 2h ( p. By (4) Lem. 4.5], we reach the desired precision v when 
hcf> > e($S/p)(v — u r ), where r is the Okutsu depth of Ftp, and v r is the Okutsu 
invariant described in section ri.il 

Therefore, the results of section [1] lead to the following routine to compute the 
different of the local extension itp /K p . 

Routine Different (char(K) = 0) 

INPUT: 

— A monic irreducible polynomial f(x) £ A[x]. 

— An OM representation tip, as in ([6]), of a prime ideal *p of B. 

1 e-*-e ei---e r) n <- J2i<j< r ( e jfj ' ' ' e rfr ~ l ) h jK e i • • ■ e j)- 

2 precision <— ev p (e), p ^— 0. 
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3 IF precision ^ THEN DO 

3.1 Apply the SFL algorithm to compute an Okutsu approximation G A[x] 
to F<p, such that > precision. 

3.2 Compute p u r +i (</>') — e-/x. (use Proposition II .5j) 

4 Dif f <- e - 1 + p. 

OUTPUT: 

— The ^P- valuation of Diff(Lrp / iTp), as the value of the variable Dif f . 

Routine Different (char(K) = p > 0) 

INPUT: 

— A monic irreducible polynomial f(x) G A[x]. 

— An OM representation t<p, as in ([5]), of a prime ideal *p of B. 

1 e «- e ei • • • e r , fi <- J2i<j< r ( e jfj ' ' ' e r/r ~ l)fy/7( e i ' ' ' e i)- 

2 p <— 0, <— </> r +i, h 4— h r +i. 

3 IF p | e THEN DO 

3.1 Compute p <— u r+ i(0') — e • p. (use Proposition II .5|) 

3.2 WHILE h < p DO 

3.2.1 Apply one iteration of the SFL algorithm to compute an Okutsu 
approximation <E> G A[x] to Fsg(x), such that /i$ > 2/i^. 

3.2.2 Compute p <— u r+ i($') — e • p. (use Proposition [L5|) 

3.2.3 h <- 

4 Dif f <- e - 1 + p. 

OUTPUT: 

— The ^-valuation of Di&(Lm/Kp), as the value of the variable Dif f . 

Let us analyze the complexity of this routine. All tasks we are interested in may 
be performed modulo p", for a sufficiently high precision v. Thus, we may assume 
that the elements of A are finite 7r-adic developments, for some it E A which is a 
local generator of p. In particular, the computation of the p-adic valuation v v has 
a negligible cost. 

Definition 2.2. An operation of A is called p-small if it involves two elements 
belonging to a fixed system of representatives of A/p. 

Let q := #A/p. A p-small operation is equivalent to O (log(g) 1+e ) word ope- 
rations, the cost of an operation in the residue field A/p. Working at precision 
v, each multiplication in A costs 0(v 1+e ) p-small operations if we assume the fast 
multiplications techniques of Schonhage-Strassen [14]. 

Denote from now on, 

n:=[L:K], e := e(qj/p), / := /(«p/p), n v := ef, <% := ^(Disc^)). 

Note that n, e, /, nqj are intrinsic invariants of L/K, but <5qj depends on the choice of 
the polynomial f(x) G A[x] used to construct this extension. In the tamely ramified 
case the routine is trivial; thus, we assume from now on that the characteristic p of 
A/p divides e. In particular, the different ideal Diff(Ly$/Kp) is divided by the e-th 
power of the maximal ideal, and n<p < v p (Disc(L<^ / K p )) < <5<p. 
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Lemma 2.3. The SFL algorithm requires O (nn<$(h/e) 1+e + n(<5<p) 1+e ) p-small 
operations in A to compute an Okutsu approximation (f> to F<$, such that h^ > h. 

Proof. Let r be the Okutsu depth of i*sp and consider the Okutsu invariants ej , fj , 
fj,j, Vj, Vj, described in section Tl. II Along the proof of Lem. 6.5], it is obtained 
an estimation of O (nn<^((h/e) 1+e + (V r +i / e) 1+e )) p-small operations for the SFL 
routine. The lemma is then a consequence of 

(7) Vr+i/e = n r + v r < 2/i r < 25<p/n<£. 

The inequality v r < p r is a consequence of e^/i > 1, for all 1 < i < r, and the 
explicit formulas for these Okutsu invariants given in section 11.11 The inequality 
Hr = n{Frp) < (%/n<p follows from Proposition II .41 and Theorem II .71 □ 

By [6] Lem. 4.21], the computation of v r+ i(<f)') is essentially equivalent to the 
computation of the ((f>i, . . . , r )-multiadic expansion of <p' . By [1] Lem. 5.4], this 
requires 0((nqj) 1+£ ) operations in A. By Remark l2.11 f2) and (J7J, this is equivalent 
to 

(8) O + K + i)/e) 1+£ ) = O ((n v ) 1+e (^/e) 1+e + (<%) 1+e ) 
p-small operations in A. 

Theorem 2.4. Assuming *}3 wildly ramified over p, the routine Different has a 

cost of O (n(<5<p) 1+e ) p-small operations in A 

Proof. Suppose first char(i^) = 0. By Lemma |2"31 step 3.1 has a complexity of 

O (nn v v p (e) 1+t + n(<%) 1+£ ) = O (n(n v ) 1+e + = O (n(6 v ) 1+i ) , 

the last equality because n<p < Srp. By (j8]), the complexity of step 3.2 is dominated 
by that of step 3.1. This ends the proof of the lemma in this case. 

If char(iir) = p > 0, the routine Different is less efficient in practice because 
we must compute iv+i (4>') after every iteration of the SFL routine. By the remarks 
following Corollarv ll.131 the number of iterations is bounded from above by log 2 (p), 
where p := p(F<$). 

By Lemma |2~51 the complexity of all SFL iterations is 

O (nn v (p/e) 1+ c + n(<5<p) 1+e ) = O . 

In the last equality we used n^pje — f p < (5<p, by Theorem 11.71 

Assume that we start the iterations with h$ — 1 (the worst possible case). If we 
compute each Okutsu approximation at the precision described in Remark 12.11 
the complexity of the computation of all tv+i ((/>') is, by flSJ): 

°f><p2V e ) 1+e + (M 1+e = O ((n<pp/ e ) 1+e + (<5<p) 1+£ log 2 p)=0 ((<%) 1+e ) ■ 

□ 

2.2. Local computation of discriminants. The discriminant of L/K is the 
product of the local discriminants, so that 

i; p (Disc(L/X)) = V v p (Disc(L v /K p )). 
*■ — 'h>Ip 

Since the local discriminant ideal Disc(Lqj / K p ) is the norm of the local different 
ideal Diff (Lep/Kp), the routine Different leads in an obvious way to a routine to 
compute v p (Disc(L/K)). This routine does not require the previous computation 
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of the discriminant Disc(/) of the polynomial f(x); actually v p (Disc(/)) may be 
deduced from the identity t> p (Disc(/)) = v p (Disc(L/ K)) + 2ind p (/), where 

(9) ind p (/) := lengthy (B ® A O p )/{A[9] ® A O p ). 

This local index ind p (/) is computed by the Montes algorithm as a by-product. 

This allows us to consider a more general routine pDiscriminant that admits 
an arbitrary monic polynomial g{x) £ A[x] as input. Before describing the routine, 
let us review some generalities on discriminants of polynomials and the role of the 
local index ind p (g) for a non-irreducible polynomial. 

Let g(x) = a n x n + • • • + a±x + ao be a polynomial of degree n, with coefficients 
in K. The discriminant of g(x) is defined as: 

DiscG?) := al n - 2 H (6, -6/) = a' 1 Res(g, g'), 

where 6\, . . . , 9 n are the n roots of g{x) in an algebraic closure, with due counting 
of multiplicities. Clearly, for arbitrary a,b £ K* , one has 

T)isc{ag(bx)) = a 2n - 2 b^- n Disc( 5 (x)). 

Thus, as regards the computation of Disc(g), we may assume from now on that 
g(x) is a monic polynomial with coefficients in A. 

Let g(x) = G\(x) ■ ■ ■ G t (x) be the factorization of g(x) into a product of monic 
irreducible polynomials in Op [x] . Let Lq be the finite extension of K p determined 
by each irreducible factor G(x), and denote by CPg the maximal ideal of the ring 
of integers of Lq- 

The discriminant has a well-known good behaviour with respect to products: 

(10) Disc(.g) = TT Diflc(Gi) TT Res(G i , Gj). 
We define the p-index of g(x) to be: 

(11) ind p (.g) := ind^d) + v P (Res(G h Gj)), 

where ind p (Gi) is, by definition, the local index ind(Gi) that was considered in 
section [TJ The reader may check that this definition coincides with that of ([9|) 
when g(x) is irreducible in A[x\. 

For each i we have w p (Disc(Gi)) = 2ind p (Gi) + v p (Disc(LGjKp)). Therefore, 
([TO]) and (UT|) show that: 

v p (Disc(g)) =2ind p (g) + V v p (Disc(L G jK p )). 

— ^is ? s* 

Now, if g(x) is a separable polynomial, the Montes algorithm computes OM repre- 
sentations of all these local factors and also the value of ind p (g) as a by-product. 
Thus, properly combined with the routine Different, we get the following routine 
to compute ^2 1<i<t v p (Disc(LGj K p )) and u p (Disc(p)). 

Routine pDiscriminant 

INPUT: 

— A monic polynomial g(x) G A[x] such that Disc(g) ^ 0. 

— A non-zero prime ideal p of A. 

1 Apply the Montes algorithm to get OM representations of the different irreducible 
factors G(x) £ O p [x] of g(x), and the value of ind p ((/). 
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2 Disc <- 0. 

3 FOR each factor G, with OM representation t^, DO 

3.1 /^/o/i •••/,. 

3.2 Call Dif f erent(<7, to) and accumulate / • w<p G (Diff (Lc/K p )) to Disc. 

4 DiscPol <- Disc + 2ind p (g). 

OUTPUT: 

— The value of Disc is the sum of the p-adic valuations of all local discriminants 
Disc(Lc/Kp), for G running on all irreducible factors of g over K p . The value of 
DiscPol is the p-valuation of Disc(g). 

Theorem 2.5. Let n := dcgg, S := v p (Disc(g)) and q ;= #A/p. The complexity 
of the routine pDiscriminant is O (n 2+e + n 1+e (l + 5) \ogq + n 1+e S 2+e ) p-small 
operations in A. 

Proof. The complexity of the Montes algorithm was analyzed in [2] , [T3] . A sharper 
estimation has been obtained in PQ, according to which, step 1 of the routine 
pDiscriminant has a cost of O (n 2+e + n 1+e (l + S) logq + n 1+e i5 2+e ) p-small oper- 
ations in A. 

By Theorem 12. 41 the complexity of step 3 may be estimated as: 

O n((5<p) 1+e ) = O (nS 1+e ) p-small operations, 

where *P runs on the wildly ramified primes lying over p. Clearly, this cost is 
dominated by that of step 1 . □ 

Corollary 2.6. If we assume p small (logq = 0(1)), then the complexity of 
pDiscriminant is O (n 2+£ + n 1+e (5 2+c ) word operations. 

In this routine the input polynomial must have Disc(g) ^ 0, otherwise the Montes 
algorithm enters into an endless loop. Nevertheless, if c\i&t(K) — 0, we can use 
upper bounds for the discriminant to design an algorithm that works for an arbitray 
monic input g(x) £ A[x\. For instance, if A = Z and p = pZ, for some prime number 
p, then Mahler's bound [5]: 

| Disc( 5 )| < n"||.g||^ 2 , := \a Q \ + ■■■ + \a n \ 

leads to v p (Disc(g)) < nv p (n) + (2n — 2)v p (\\g\\ 00 ). The reader may derive similar 
upper bounds, w p (Disc(f/)) < N(g), in the general case. Along the flow of Montes 
algorithm, partial p-indices (taking positive integer values) are accumulated to a 
variable Index, whose final output value is ind p (<7). If Disc(g) = 0, we run into an 
endless loop that strictly increases Index at each iteration. Thus, we may introduce 
a control instruction that allows the next iteration while the value of Index is less 
than N(g), but it breaks the loop and outputs i>p(Disc(g)) = oo otherwise. 

3. Computation of the p-adic valuation of resultants 
Consider two polynomials 

f(x) = a n x n + • • • + a\X + a , g(x) — b m x rn + • • • + b\x + b £ K[x], 
of degree n, m, respectively. For any a, 6, c £ K* we have 

Res(af(bx),cg(bx)) = a m c n b nm Res(f(x), g(x)). 
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Hence, as regards the computation of the resultant, we may suppose that f(x),g(x) 
are monic and have coeficients in A. 

In [SI Sec. 4.1] the p-adic valuation of the resultant of two monic polynomials 
in A[x] is computed by applying a kind of non-optimized version of the Montes 
algorithm. In this section we reformulate those ideas into a concrete algorithm 
based on the optimized Montes algorithm as described in [5]. 

3.1. Partial exponents of resultants. In this section we introduce certain pos- 
itive integers Res t ^(/, g), whose accumulation is equal to u p (Res (/,<?)). These 
partial u p -values of resultants are computed in terms of combinatorial data of ade- 
quate Newton polygons of the polynomials f(x), g(x). 

Let us describe in some detail the relevant Newton polygon routine. Along the 
algorithm that computes w p (Res(/, <?)), we shall construct types of order i — 1, 

t = Oo; (0i, Ai, V'i); • • • ; (<fo_i,Ai-i>V ; i-i))> 

dividing f(x). All polynomials cf>i, . . . , (f>i—i will have coefficients in A, and the 
type will always be loaded with three level i invariants, <f>i, csi, uji. The invariant 
4>i is a representative of t; that is, 4>i(x) € A[x] is a monic polynomial of degree 
rrii := ei_i/j_imj_i, such that i?j„i(</>i) ~ i/'i-i- The choice of a representative of 
t determines a Newton polygon operator Ni, with respect to <pi and the MacLane 
valuation Vi supported by t (see section [O]) . 

The invariant cs^ is a cutting slope; it is a positive integer h telling us that we 
must compute only the subpolygon N^(f) formed by the sides of ATj(/) of slope less 
than — h. By a variation of Lemma 1 1.11 we are able to compute a priori the length 
of this subpolygon, and this is the role of the third invariant: uJi = £(N^ Si (/)). 
Since we always know a priori the length of the Newton polygons we are interested 
in, we may use the following Newton routine. 

Routine Newton(t, uj, g) 

INPUT: 

— A type t of order i — 1 > and a representative <fi of t. 

— A positive integer to. 

— A non-zero polynomial g{x) G i4T[a;]. 

Compute the first uj + 1 coefficients ciq(x), . . . , a LU (x) of the canonical (^-expansion 
of g(x) and the Newton polygon N of the set of points (s, Vi(a s (/> s )), for < s < uj. 

OUTPUT: 

— A list [Si, . . . , S t ] of the sides of N. 

A side S of N is a segment of negative slope of the Euclidean plane, whose 
end points have non-negative integer coordinates. The length, E = £(S), and 
height, H , of S are the lengths of its projections to the horizontal and vertical axes, 
respectively. Following the convention of [51 Sec. 1.1], we admit sides S of slope 
— oo; we may think that the left end point of such an S is (0, oo), and the right end 
point is (s, u), with s € Z>o and u G Z>o. The length of such a side is E = s and 
the height is H = oo (see Figure [3]). 

If the left end point of a Newton polygon N has a positive abscissa, we consider 
(formally) that the side of slope — oo determined by this point is also one of the 
sides of TV. According to this formalism, if N = Si + • • ■ + St are all sides, finite 
and infinite, of N, we have always £(N) = £(Si) H h £(S t ). 
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Figure 3. Length and height of a side of negative slope 




Definition 3.1. Let S,S' be two sides of negative slope, with lengths E 7 E' and 
heights H, H' , respectively. Let h be a non-negative integer. We define 

Res(S, S') := min{EH' , E'H}, Res h (S, S') := Res(S, S') - hEE'. 

Let t = (tpo; (<pi, Ai, -0i); • ■ ■ ; (<fri-i, Xi-i, V'i-i)) ^ e a tyP e °f order i — 1 and Zet 4>% 
be a representative oft. Let f{x),g{x) € A[x] be two monic polynomials. We define 

Res? a- (/> ff) := /o ■ ■ • /i-i Res k (S, T) 

= fo • • ■ /i-1 

As the notation suggests, Res t ^ (/, gr) depends on the choice of the representative 
4>i . The value of Res£ ^ (/, g) is equal to infinity if and only if both Newton polygons 
Ni(f), Ni(g) have a side of slope -co. This happens if and only if both polynomials 
/, g are divisible by 4>i m A[x\. 

If Res^ i (/,<?) is finite, then it takes a non- negative integer value. It vanishes if 
and only if either Ni(f) or Ni(g) have no side (neither finite nor infinite) of slope 
less than —h; that is, if and only if either N^(f) or N^(g) have length zero. 

3.2. The algorithm. Roughly speaking, the algorithm consists of a simultaneous 
application of the Montes algorithm to f(x) and g(x), and the accumulation of all 
values Resjj s ^. (/, g), for all the types t considered along the flow of the algorithm, 
such that t divides both polynomials f(x) and g{x). 

The Montes algorithm is fully described in [5], in terms of the theoretical back- 
ground developed in [6]. For a short review, we address the reader to [TO] , or [TJ 
Sec. 4]. 

We present a detailed description of the routine pResultant, directly in pseu- 
docode. We denote by <f)\, A*, Vf, etc. the data at the i-th level of a type t. 

Routine pResultant 

INPUT: 

— Two monic polynomials f(x),g(x) G A[x], with Res(/,<?) ^ 0, deg/ < degg. 

— A non-zero prime ideal p of A. 



E 



Res(5,T) 



SeNNf),T£NNg) 



he(N?(f))£(N t h (g)) 



1 Set Res value <— 0. 
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2 Factorize f(y) (mod p) = Y\ ip p(y) a " ( ' over (A/p)[y\. 

3 FOR every irreducible factor cp DO 

4 Set b <— oid v (g (mod p)). IF b = THEN continue to the next factor cp. 

5 Take a monic 4>(x) £ such that 0(y) (mod p) = <p(y). 

6 Create a type t of order zero with: ip^ ^— cp, <fi\ <— <— a v , a* <— b, and 
cs* •<— 0. Initialize the list Stack <— [t]. 

WHILE Stack is non-empty DO 

1 Extract (and delete) the last type to from Stack. Let i — 1 be its order. 

2 Call Newton(to, u>i, f(x)) and Newton(to, cti, g(x)). 

3 Compute Res^ 1 ^. (/, g) and add this number to Resvalue. 

4 FOR every finite side S of A^ cs< (/) DO 

5 Set A* <- slope of 5. IF A* is not a slope of N™ 1 (g) THEN continue 
to the next side S. 

6 Compute Ri(f),Ri(g) G F 4 [y]. Factorize i?,(/) over F<[y]. 

7 FOR every monic irreducible factor ip of Ri{f) DO 

8 Set 6 <- ord^ R l (g). IF 6 = THEN continue to the next factor ip. 

9 IF oji > 1 THEN make a copy t of the type to and extend it to an 

order i type by setting: ip\ <— ip, <— ord^, Ri(f), <— b. 

Compute a representative <p £ A[x] of t. 

ELSE apply one iteration of SFL to compute with h$ > Ih^t . 

10 IF deg 4> = deg$ THEN set 

4>\ <- 4>, cs* <r- \X\\, u\ <- <4+i, a\ <- a\ +1 , 
and delete all data in the (i + l)-th level of t 
ELSE set <p\ +1 <- cp, cs* +1 «- 0. 

11 Add t to the Stack. 
END WHILE 

OUTPUT: 

— The p-adic valuation of Res(/, g), as the value of the variable Resvalue. 

Remark 3.2. (1) We mentioned already that oji = i(N^ Si (f)). Thus, the compu- 
tation (and storing) of this invariant saves operations in the Newton routine because 
we know a piori how many coefficients of the (^-development must be computed. 
For the same reason, we compute and store in all types a completely analogous 
invariant on telling a priori the value of £(N^ Si (g)). 

(2) In step 9, the property u)i = 1 implies that to is already /-complete (Defini- 
tion [TT^]), and <p\ is an Okutsu approximation to one of the irreducible factors (say) 
F(x) of f(x) over O v . In this case, the SFL routine may be applied, and one single 
iteration of its main loop leads to an Okutsu approximation (p with > [8]. 
On the other hand, the standard construction of a representative of the type t leads 
to an Okutsu approximation <p with > h^t. Thus, the use of SFL accelerates the 
process of getting <fi sufficiently close to F, to have (to; (<fi, A, ip)) \ g(x) in a future 
iteration of the WHILE loop. 
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This acceleration leads to a significant improvement of the routine in cases where 
the polynomials f(x) and g(x) have p-adic irreducible factors which are very close 
one to each other. 

(3) Step 10 takes care of the optimization. If deg<j) = deg^*, then the future 
information provided by the Newton polygon JV i T 1 (/), with respect to the pair 
(t,</>), is equivalent to the information provided by N^ Si (f), with respect to the 
pair (to,</>) [HI Sec. 3.2]. The latter option is more efficient because we work at a 
lower order. 

Theorem 3.3. The routine pResultant terminates and its output value is in- 
deed u p (Res(/, g)). It requires O (n 2+e + n 1+e (l + 5) \ogq + n 1+e i5 2+e ) p-small ope- 
rations in A, where n = max{deg<?,deg/i}, S = i>p(Res(g, h)) and q = fj=A/p. 

Proof. Let us call basic pResultant the algorithm obtained by eliminating the 
optimization procedure of step 10; that is, by replacing this step by: 

10' Set (t>\ +l <- (f>, cs* +1 <- 0. 

By O Thm. 4.10], the basic pResultant algorithm terminates and the final 
output value of the variable Resvalue is u p (Res(/, <?)). 

On the other hand, pResultant and basic pResultant are equivalent algo- 
rithms, in the sense that they have the same number of iterations of the WHILE 
loop and they accumulate the same value to the variable Resvalue at each ite- 
ration. This is a consequence of [5J Props. 3.4+3.5] and the rest of the arguments 
of [3 Sees. 3.2+3.3], where a completely analogous situation was discussed when 
we compared basic and optimized versions of the Montes algorithm, both leading 
to the computation of ind p (/) by the accumulation of partial indices. 

The complexity analysis is obtained by completely analogous arguments to those 
used in [1] to estimate the complexity of the Montes algorithm. □ 

Corollary 3.4. If p is small, then pResultant requires O (n 2+e + n 1+e <5 2+e ) word 
operations. 

As it was the case for the routine pDiscriminant, if char(K) = 0, the assumption 
Res(/, g) may be omitted by considering an upper bound for i?p(Res(/, g)). For 
instance, for A — Z, n — degf <m = degg, we have [3J Thm. 6.23], 

Res(f, g ) < \\m\\9E = (kol 2 + • • ■ + KI 2 ) m/2 (M 2 + • • ■ + IM 2 )" /2 . 

We leave to the reader the derivation of similar upper bounds, u p (Res(/, <?)) < 
N(f, g), in the general case. If Res(/, g) — 0, the WHILE loop of pResultant is an 
endless loop. However, it increases Resultant by a positive integer value at each 
iteration. Thus, we may introduce a control instruction allowing the next iteration 
as long as the value of Resultant is less than N(f, g), and breaking the loop with 
the output u p (Res(/, g)) = oo otherwise. 

4. Numerical tests 

We have implemented the routines pDiscriminant and pResultant in the case 
A = Z. They are included in the Magma package pFactors.m, which may be 
downloaded from http : //montespro j ect . blogspot . com. 

We present in this section some numerical tests for several polynomials selected 
from the families of test polynomials in the appendix of [8]. All tests have been 
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done in a Linux server, with two Intel Quad Core processors, running at 3.0 Ghz, 
with 32Gb of RAM memory. Times are expressed in seconds. 

We compare running times of our routines (abbreviated as pDis, pRes in the 
tables), with the naive routine that computes first Disc(g) or Res(<?, h), and then its 
p- valuation. The numerical results show that the pDiscriminant and pResultant 
routines are much more efficient than the corresponding naive routines, when the 
degree of the polynomials and/or the size of the coefficients grow. 

4.1. Numerical tests of pDiscriminant. 

Example 1. Let p be a prime number and n a positive integer. Consider the 
polynomial of degree n: 

f(x) = {x+p + p 2 + ---+p 20 ) n +p 20n+1 . 

This polynomial is irreducible over Z, p [x] and it has Okutsu depth equal to one. It 
has u p (Disc(/)) = (20n + l)(n - 1) + nv p (n) App.]. 



p 


deg/ 


«„(Disc(/)) 


pDisc 


naive 


2 


20 


7659 


0.01 


0.00 


2 


50 


49099 


0.02 


0.83 


2 


100 


198299 


0.07 


26.60 


2 


150 


447299 


0.13 


190.57 


59 


20 


7619 


0.02 


0.00 


59 


50 


49049 


0.06 


66.14 


59 


100 


198099 


0.22 


1887.55 


59 


150 


447149 


0.59 


12225.08 


59 


500 


4990499 


12.80 


> 24 hours 



Example 2. Let us test the influence of the existence of many p-adic irreducible 
factors. Let p > 5 be a prime number, and m a positive integer, m < p/2. Let 
go(x) = x w + 2p n . Consider the polynomial of degree 10m: 

g(x) = g (x)g (x + 2) ■ ■ ■ g (x + 2(m - 1)) + 2p 110m . 

This polynomial has m irreducible factors of degree 10 over 1 p [x\, all of them with 
Okutsu depth equal to one. It has v p (Disc(g)) — 99m 8., App.]. 



p 


m 


deg g 


v p (Disc(g)) 


pDisc 


naive 


7 


3 


30 


297 


0.01 


0.22 


23 


5 


50 


495 


0.03 


8.30 


23 


10 


100 


990 


0.04 


243.76 


43 


5 


50 


495 


0.04 


12.92 


43 


10 


100 


990 


0.04 


378.17 


43 


15 


150 


1485 


0.04 


2598.82 


101 


50 


500 


4950 


0.27 


> 24 hours 



Example 3. Finally, let us test the influence of the growth of the Okutsu depth in 
both routines. Let p > 3 be a prime number. Consider the following polynomials: 
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E 1 (x) 


= x 2 + p 


E 2 {x) 


= E^xf + (p-l)p 3 x 


E 3 (x) 


= E 2 (xf+p u 


E 4 (x) 


= E 3 {xj 3 +p 29 xE 2 (x) 


E 5 (x) 


= Ei{xf + (p- l)p i ' 1 xE l {x)E 3 {x) 2 


E 6 (x) 


= E 5 (x) 2 +p 88 xE 3 (x)E i (x) 


E 7 {x) 


= E 6 (x) 3 +p 295 E 2 (x)E 4 (x)E 5 (x) 


Ea(x) 


= E 7 (x) 2 + (p- l)p 632 xE 1 (x)E 2 {x) 2 E 3 (x) 2 E 6 {x) 



These polynomials are all irreducible over Z p [x] and the Okutsu depth of Ej is 
equal to j, for all 1 < j < 8 |U App.]. 



p 


j 


deg Ej 


v p (Disc(Ejj) 


pDisc 


naive 


5 


5 


72 


4671 


0.01 


0.09 


5 


6 


144 


18899 


0.02 


1.66 


5 


7 


432 


171383 


0.15 


204.54 


5 


8 


864 


686825 


0.72 


4565.77 


61 


5 


72 


4671 


0.01 


0.39 


61 


6 


144 


18899 


0.02 


9.63 


61 


7 


432 


171383 


0.22 


1580.66 


61 


8 


864 


686825 


1.25 


35623.63 



4.2. Numerical tests of pResultant. 

Example 4. Let p > 5 be a prime number, m a positive integer, m < p/2, and 
g{x) the polynomial of Example 2. Let f(x) be the polynomial of Example 1 of 
degree n = 10m = deg g. 



p 


deg / = deg g 


i> p (Res(/,ff)) 


pRes 


naive 


7 


30 


300 


0.01 


0.94 


11 


50 


500 


0.00 


19.58 


23 


100 


1000 


0.01 


1060.03 


31 


150 


1500 


0.01 


8433.39 


43 


200 


2000 


0.03 


38430.25 


101 


500 


5000 


0.27 


> 24 hours 



Example 5. Let p > 3 be a prime number, and consider the polynomials Ej(x), 
1 < j < 8, of Example 3. 
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if 


% 


7 




pRes 


naive 


5 


5 


6 


9557 


0.04 


0.23 


5 


5 


7 


28671 


0.14 


2.70 


5 


5 


8 


57342 


0.61 


13.76 


5 


6 


7 


57343 


0.15 


13.11 


5 


6 


8 


114686 


0.67 


66.12 


5 


7 


8 


344059 


0.83 


885.89 


101 


5 


6 


9557 


0.05 


1.95 


101 


5 


7 


28671 


0.22 


28.11 


101 


5 


8 


57342 


1.12 


150.30 


101 


6 


7 


57343 


0.25 


141.52 


101 


6 


8 


114686 


1.31 


746.78 


101 


7 


8 


344059 


1.60 


9335.33 
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